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Abstract 

We analyze the effect of a colored non Gaussian noise on a model of a random walker moving along 
a ratchet potential. Such a model was motivated by the transport properties of motor proteins, like 
kinesin and myosin. Previous studies have been realized assuming white noises. However, for real 
situations, in general we could expect that those noises be correlated and non Gaussian. Among 
other aspects, in addition to a maximum in the current as the noise intensity is varied, we have also 
found another optimal value of the current when departing from Gaussian behavior. We show the 
relevant effects that arise when departing from Gaussian behavior, particularly related to current's 
enhancement, and discuss its relevance for both biological and technological situations. 

PACS numbers: 05.40.Jc, 82.20.Uv,05.40.-a 
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I. INTRODUCTION 



Noise induced transport by Brownian motors or "ratchets" has attracted the attention 
of an increasing number of researchers due to its biological interest as well as its potential 
technological applications. Since the pioneering works, besides the built-in ratchet-like bias 
and correlated fluctuations, several different aspects have been studied, such as tilting 1, 2| 
and pulsating potentials, velocity inversions 1, 4], etc. There are relevant reviews 5, 6| 
indicating the biological and/or technological motivations for the study of ratchets as well 
as showing the state of the art. 

Among other aspects, ratchets has been used to explain the unidirectional transport of 
molecular motors within a biological realm [s], [g]. Among the different motor proteins, ki- 
nesin has attracted considerable attention motivated by experimental results in which the 
dynamical details of its motion can be measured Kinesin is a protein with two 

heads that performs a walk along the microtubule inside cells. Motivated by these experi- 
mental results, several researchers have introduced diverse models in order to understand the 



particular form of walking of kinesin 



10|. Usually those models consider a walker moving 



along an asymmetric ratchet potential mediated by noise. This walker has two feet that 
are represented by two particles nonlinearly coupled through a bistable potential. Among 



11 



12| including all the above indicated 



those models, there is a recent one introduced in 
ingredients, where the walker moves along a track formed by an asymmetric potential, being 
subjected to two independent white noise sources acting on each of the two particles and to 
a common external harmonic force. It was observed that the current J as function of noise 
intensity presents a maximum, conforming another example of the constructive role of 
noise. 

Recent studies on the effect of a non Gaussian noise on several noise induced phenom- 
ena, have shown the existence of strong effects on the system's response. In those studies, 
opposing to the most usual cases that only consider Gaussian white noises, such an effect 
was analyzed in stochastic resonance, noise induced transitions, "standard" ratchets, etc 
This form of noise was motivated by the nonextensive statistical me- 
19|. In this work we want to analyze the effect of this form of colored and non 



n 



chanics [l8|, ^ 

Gaussian noise over the kinesin ratchet model introduced in Refs. 
pose, we used a mean field approximation and exploited a recently developed technique [l7| 



\\\ . For this pur- 
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Through the variation of a parameter g, this form of noise offers the possibihty of analyzing 
the departure from Gaussian behavior (corresponding to g = 1). Since subtle change of 
environment conditions can produce drastic changes in biological process, such effects could 
be very relevant. It is worth here indicating that a related mean field approximation was 
introduced in [2^. 

The organization of the paper is as follows. In the next Section we show the mean-field 
approximation and how reliable it is, based on the good agreement with known results for 
the Gaussian case. The third Section presents the statistical properties of the non Gaus- 
sian noise, and the form of the current that results after applying the effective Markovian 
approximation in the present problem. In the following Section we present the results when 
we depart from the Gaussian behavior, that is varying the parameter q. Finally, in the last 
Section we present some conclusions. 



II. THE MEAN-FIELD APPROXIMATION AND ITS ACCURACY: GAUSSIAN 
CASE 



The stochastic dimensionless differential equations for the two particles in the overdamped 
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regime, whose coordinates are indicated by x and y, are 

X = -d^V{x) - -y) + \l^wiiit) + A sin fit, 

y = -dyV{y) - dyVbix - y) + ./2I\,^2it) + A sin fit, 

where V{x) is the dimensionless ratchet potential 

V{x) = C + Ub^ [sin(27r(a; - a;o)) - 0.25 sin(47r(a; - a;o))] 



(1) 



The constant C = UFi[sin{2TTXo) — 0.25 sin(47ra;o)] was chosen in such a way that V{0) = 0. 
The constant Xq is introduced in order to center the minima of the periodic potential on 
integer values. The dimensionless amplitude of the ratchet potential is indicated by Ur. 

The dimensionless bistable potential Vb{x — y), that represents the nonlinear coupling 
between the two particles, is given by 



Vbix -y) = Ub 



(x - yY _ 2 - yf 



(2) 



/4 P 

where Uh is the dimensionless amplitude of this bistable potential and 21 is the distance 
between the two minima. 
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In the original model, it is assumed that ^i(t) and ^2(^) are Gaussian white noises with zero 
mean and correlation {C,i{t)^i{t')) = 2Duj6ij6{t — t'), with the intensity of the statistically 
independent noises. In our model, as indicated at the introduction, we assume they are non 
Gaussian colored noises, with characteristics that we will briefly indicated later. Clearly, 
the model also considers an external harmonic force. 

In [12I this model was analytically solved. However, in order to consider the non Gaussian 
colored noise case, we need to introduce a different approach. For that purpose, we consider 
a mean field approximation (MFA) 2l(|. For the present case, such a MFA consist in the 
following approximation for Vb{x — y) 



Vb{x -M) = Ub 



X - Mf (x - Mf 



l + - ^^-2 



(3) 



/4 P 

with M =< y >, and clearly we have < x >=< y >. After applying the MFA, the equations 
for both variables x and y results to be of the same form. Hence, we can reduce the problem 
to a single equation describing the model system 



X = M, t) + pD^i{t), (4) 

where Veff{x, M, t) = V{x) + Vb{x - M) - xA sin Qt. 

In order to test the MFA, we start analyzing the Gaussian white case, and show that it 
gives similar qualitative results than the numerical simulations and the analytical solution 
for the original model 11, 1^. Hence, in Eq. (jlj), we start assuming that ^{t) is a Gaussian 



white noise with the same behavior than the ^j(t)'s. 

Considering an adiabatical approximation in order to decouple the external harmonic 
force from the rest, the corresponding Fokker-Planck equation can be solved assuming pe- 
riodic boundary conditions. Even though Veff{x,M,t) is non periodic [due to Vb{x — M) 
being non periodic], we can still assume such a periodic behavior since both feet are never 
too far away from each other, the strong slope of the attractive potential preventing it to 
occur. In Fig. 1 it is shown Vb{x) vs. x with and without the indicated approximation. In 
Fig. 2 we show Ve//(x) vs. x, again with and without the approximation. It is apparent 
that the curves looks quite alike. 

Considering the adiabatical approximation, the corresponding stationary solution of the 
Fokker-Planck equation is 22] 
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FIG. 1: The bistable potential Vb vs. x. The continuous line corresponds to the exact form of the 
potential, while the dashed line is for the approximate one. The parameters are Ub = 0.2512 and 
I = 1. 



with 



H(x) 



1 



x+L 



a normalization constant, and L = 2tt the period. Hence, numerically solving 

rL/2 

we can obtain the current J{t) [22] as 



l-L/2 

/ dxxP'\x,M) = M =< X > 



1 



,{V,ff{L)-V,f!m/D^ 



and the net current is obtained as 



with T = 



(6) 



(7) 



(8) 



(9) 



In order to check these results in Fig. 3 we depict J vs. Dw The curve presents a 
maximum for an "optimal" value of the noise intensity, which is qualitatively similar to the 



one obtained in 



11 



12j. With the support offered by this agreement, we are now in position 



to use this same approach for the non Gaussian colored noise case. 




X 



FIG. 2: Effective potential Vgff vs. x. The continuous line corresponds to the exact form of the 
potential, while the dashed line is for the approximate one. The parameters are A = 1, I = 1, 
Ur = 0.16, Ub = 0.2512 and t = 0. 

III. NON GAUSSIAN NOISE CASE 



A. Statistical properties of the non Gaussian noise 



We consider now the non Gaussian colored noise case as in 



13|. Hence, in Eq. ^{t) 



is a noise with a dynamics described by the following Langevin equation 

m = • 

with ri{t) a white Gaussian noise, and 



(10) 



K(0 



A, 



In 



(11) 



For q = 1, the process ^ coincides with the Ornstein-Uhlenbeck (OU) one (with a correlation 
time equal to r), while for g 7^ 1 it departs from the Gaussian behavior. For g < 1 the 
stationary probability distribution (spd) has a bounded support, with a cut-off given by 
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FIG. 3: Net current J for a gaussian white noise, as function of D^. The values of parameters are 
A = l,l = l,Ur = 0.16 and Ub = 0.2512. 



Il^ll = CO = [{1 — q)T/ {2D^)] 2 J with a form given by 



1-9 



(12) 



for ll^ll < uj and zero for ||,^|| > uj {Zg is a normalization constant). Within the range 
1 < g < 3, the spd is given by 



1 + 



2A„ 



1-9 



(13) 



for — oo < ^ < oo, and decays as a power law (that is, slower than a Gaussian distribution) 
for ^ ^ oo. Finally, for g > 3, this distribution can not be normalized. Note that keeping 
constant, when increasing q the dispersion of the distribution also increases. In 13|] the 



second moment of the distribution was obtained. This moment is related to the intensity of 
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the non Gaussian noise, and is given by 



D. 



ng 



--< r >= 



2D„ 



(14) 



r(5-3g)^ 

that diverges for g > 5/3. For Tng, the correlation time of the process C,{t), that was defined 
in detail in [13], it was not possible to find an analytical expression. However, it is known 
Q that fo.;: 5/3. dive, ges as (5 — 3g) ^. In [l3| Tng vs. q was numerically calculated 
for the range 0.5 < g < 5/3. The following analytical approximation, accurate within the 
indicated range, was found 

2r 



r, 



ng 



(5 - 3g) 



[l + 4(g-l) 



(15) 



B. Current in the non Gaussian case 

The Fokker-Planck equation for the present case, that is the system described by the 
Langevin equation indicated by Eq. (jlj) with ^(t) the non Gaussian noise, has a known 
structure. After applying the effective Markovian approximation as exploited in 13. [isl. fl?!. 
such an equation has the form 

1 



where 



with 



and 



with 



dtP{x,t) = d^[A{x)P{x,t)] + ^dl[B{x)P{x,t)], 



A{x) 



V'fr 



Ai{x) 



Ai{x) + A2{x) 

l-(r/2iPj(g-l)\//|/ 
l + ir/2D^)iq-l)V:jf 



A,ix)=TV:'ff[l + iT/2D^){q~l)V, 



Sf^' 



B(x) 



[l + {T/2D^){q-l)V:' 1212 



effi 



B,{x) + B2ix) 



(16) 



(17) 



(18) 



B,ix) = rV:ff \l + ir/2D^)iq-l)V, 



/2 

eff 



Boix) 



l-(r/2Dj(g-l)K'J/ 



The prima indicates a derivative respect to x. 



The corresponding stationary solution of the indicated Fokker-Planck equation is 22 1 



with a normahzation constant, and 



where S{x) = \fB and R{x) = —A — S'. Then, using Eqs. (JTj), as and we can 
obtain the current J. 

IV. RESULTS 

In what follows we present several results for the case of submitting the system to a non 
Gaussian noise, mainly restricting ourselves to the range 0.5 < g < 1.5, which is the range 
where our evaluation is, in principle, valid. In all cases, except we indicated something 
different, we adopted r„g = O.Olvr, A = l,l = l,n = 0.5, Ur = 0.16 and Ub = 0.2512. 

In Fig. 4 we show the dependence of J on q, for different values of Dng- Looking at the left 
part of the picture, we identify the curves, from top to bottom, with the values Dng = 0.02, 
0.038, 0.06, 0.07, 0.085, 0.1 and 0.13. The entanglement between the dependence on q and 
Dng is apparent. For low values of Dng {Dng ~ 0.02) the behavior of J corresponds to an 
initial plateau for g < 1 and a fast (step like) decrease for g > 1. When Dng increases 
{Dng ^ 0.06), J{q) adopts a bell's shape, with a maximum that shifts to larger values of q 
for increasing Dng- It is worth noting that for large values of Dng, the relevant values for 
the current will only occur when q takes values within the indicated bell-like region 
As in previous studies on the effect of non Gaussian noise on other systems 
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16l . Il7l | we conclude that the observed change of behavior is associated with the process ^(t), 
that is with the departure from Gaussian behavior (g = 1) and the associated changes in 
Pq{^). As indicated before, for g < 1 there is a cut-off that prevents that ^(t) reaches large 
values, while for g > 1 ^{t) can reach very large values. For small values of Dng it is expected 
that the current J(g) will not be very sensitive to such a cut-off, as the probability for ^{t) 
to reach those cut-off values is very small. In opposition, for large values of Dng we expect 
that J(g) becomes highly sensitive to the cut-off. It seems to be the case as J varies very 
slowly with g (for g < 1) when Dng is small, while for large values of Dng just the opposite 
occurs. 
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FIG. 4: Net current J vs. q, for the range .5 < g < 1.5. Looking at the left part of picture, we 
identify from top to bottom Dng = 0.02, 0.038, 0.06, 0.07, 0.085, 0.1 and 0.13. Other parameters 
are: r„g = O.OItt, A = l,l = l,Ur = 0.16 and Ub = 0.2512. 

However, for g > 1, it is the tail of the spd that mainly contributes to build the current 



17(1 . As indicated above, for g < 1 (that is, for the case without tail), and large values of 
Dng, J drastically decreases when q decreases. Let us select one of the curves for J with 
a bell shape, and look into the spd for values of q below, at, and above of the maximum 
of that curve. Such analysis is shown in Fig. 5. We have chosen Dng = 0.085, q = 0.94 
(below), q = 1.2 (at the maximum) and q = 1.4 (above). We observe that for the case 
of the maximum, the spd is symmetric with two peaks at both ends of the space period. 
However, when departing from such "optimal" condition, the spd looses its symmetry. In 
one hand, when q is at the left of the maximum (g = 0.94) a new peak emerges at the 
central zone of the period. On the other hand, when q is to right of the maximum, there 
are two peaks that arise between the previously indicated ones. In both cases, the prize 
of these additional peaks is that the spd decreases at both ends, indicating a reduction of 
J. The same behavior was also observed in all cases where J{q) adopts a bell's shape. We 
also compare the spd behavior for the case when J(g) has a step-like shape, and note that 
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FIG. 5: Stationary probability distribution P for q below, at, and after of maximum J's value, 
with Dng = 0.085 and: dot line for q = 0.94, solid line for q = 1.2 and dash line for q = 1.4. Other 
parameters are: r„g = O.OItt, A = 1, I = 1, Ur = 0.16 and Uf, = 0.2512. 

for values of q within the plateau, the spd is like the one corresponding to the case of the 
maximum of J{q) (Fig. 3, for q — 1.2), while for values of q beyond the step, the spd is 
analogous to the one corresponding to a value of q located at the right of the maximum 
(g = 1.4). It seems that the current reaches a maximum when one feet follows the other, 
and their separation is kept constant, and equal to a period. We have also analyzed the 
Gaussian case {q — 1) for different values of Dng, with results that essentially reproduces 
those of Fig. 5. In short, for small Dng the curves are analogous to those obtained for large 
q (having four peaks), while for large Dng the shapes resembles the ones for small q (having 
three peaks). 

Figure 6 depicts J{t) -the current before its time averaging- as function of time, for a 
fixed Dng and different values of q. It is apparent that the current adopts both, positive 
and negative values along the period of the external forcing. We observe that for q < q^ax 
(where J{qmax) reaches its largest values) there is a balance between positive and negative 
current's values yielding, in average, a small net current. When q increases, the negative 
values of the current reduces, the indicated balance is also reduced, and the net current 
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FIG. 6: J{t), the current before the time averaging, vs. t, for Dng = 0.085 and different values of 
g's. Adopting as reference the left part of the figure, from top to bottom we have q = 0.3, 1., 1.2, 
1.4 and 1.55. Other parameters are: T„,g = O.Olvr, ^ = 1, / = 1, ^7^ = 0.16 and Ub = 0.2512. 

increases. For even larger values of g (g > qmax), the positive value of the current is reduced, 
while the negative one tend to zero, and the net current decreases. 

We finally analyze the system's behavior with r„g. Firstly, we observed that inside the 
studied range, J varies almost linearly with r„g. Since the corresponding slope change its 
sign with q and Dng, we choose to look into the behavior of dJ/dTng. Hence, in Fig. 7 we 
show dJ/dTng vs. q for different values of D^g, while in Fig. 8 we show dJ/dTng vs. D^g for 
different values of q. The sensitivity of dJ/drng against variations of q and Dng is apparent, 
and the changes in the slope's sign are clearly seen. 



V. CONCLUSIONS 



In line with some recent work Il3l 



14, a 



16 



171 , we have here analyzed the effect of 



11 



12|. 



a colored and non Gaussian noise over the kinesin ratchet model introduced in 
As was discussed in 13|, there are strong experimental evidences that noise sources within 
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FIG. 7: dJ/dTng vs. q for different Dng- Considering the left part of the figure, we identify from 
top to bottom, Dng = 0.1, 0.05, 0.038, 0.03 and 0.02. Other parameters are: r„g = O.Olvr, A = 1, 
l = l^Ur = 0.16 and Ub = 0.2512. 

biological systems can not be in general Gaussian This fact gives strong support to the 
study of the effect of non Gaussian noises within a biological motivated context. 

In order to make an almost analytical treatment, we have here exploited a mean- field 



11 



12 



like approach, that shows a nice agreement when it is tested against known results 
for the Gaussian case {q = 1). Hence, we can trust the results we have obtained when 
varying the parameter q. For the adopted form of noise, such variation offers the possibility 
of analyzing the departure {q ^ 1) from the Gaussian behavior (g = 1). 

When analyzing the dependence of the current J on Dng and q, we have confirmed the 
existence, for a fixed value of q, of a maximum or "optimal" current as function of Dng- 
What is new, is that for fixed Dng, there is also an optimal value of q yielding a maximum 
value of the current. More, in general such a value corresponds to a non Gaussian situation 
(g 7^ 1). Such a results could be understood from the behavior of the spd for different values 
of q, as well as from the current's time dependence (that is, before the time averaging). 
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FIG. 8: dJ/dTng vs. Dng for different values of q: squares q=0.52, circles (7=0.85, triangles g=1.015, 
diamonds (7=1.125, stars (7=1.235 and crosses g=1.4. Other parameters are: T„g = O.Olvr, ^4 = 1, 
l = l^Ur = 0.16 and Ub = 0.2512. 



Regarding the dependence of J on r„g, we have seen that it is linear. However, analyzing 
its slope {dJ/ dTng) it becomes apparent the possibilities of even changing the sign of such a 
slope varying both, q for fixed Dng, or Dng for fixed q. 

The indicated results show us the richness of behavior that we can found when departing 
from the Gaussian situation. It is apparent that those results could be of relevance not only 



for biological studies, but for technological applications as wel 



. This point was discussed in 



171]. From a technological point 



general in [l^, and particularly for the case of ratchets in 
of view, such a noise source offers a variety of forms of controlling and/or optimizing the 

transport process. ^ ^ 

Clearly, the exploitation of an approach like the one used in [l2l], based in a coordinate 



separation, could offer a more transparent and better description of this system's physics. 
In addition a numerical validation of the whole approach is required. Also, another aspect 
that requires further analysis is the case when 21 is non-commensurate with the unit cell of 



14 



the ratchet potentiaL All these aspects will be the subject of forthcoming work 



2^ 
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